Low-momentum nucleon-nucleon interaction and shell-model calculations 
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We discuss the use of the low-momentum nucleon-nucleon (NN) interaction V[ ow _k in the deriva- 
tion of the shell-model effective interaction and emphasize its practical value as an alternative to 
the Brueckner G-matrix method. We present some selected results of our current study of exotic 
nuclei around closed shells, which have been obtained starting from the CD-Bonn potential. We also 
show some results of calculations performed with different phase-shift equivalent NN potentials, and 
discuss the effect of changes in the cutoff momentum which defines the V[ ow _k potential. 

PACS numbers: 



I. INTRODUCTION 

In recent years, shell-model calculations employ- 
ing realistic effective interactions derived from modern 
nucleon-nucleon (NN) potentials have produced results 
in remarkably good agreement with experimental data 
for a number of nuclei in various mass regions. A main 
difficulty encountered in this kind of calculations is the 
short-range repulsion contained in the bare NN potential 
Vnn, which prevents its direct use in nuclear structure 
calculations. This difficulty is usually overcome by re- 
sorting to the well-known Brueckner G-matrix method. 

Recently, a new approach has been proposed [J Q 
which consists in deriving from Vnn a low-momentum 
potential Vow-k defined within a given cutoff momen- 
tum A. This is a smooth potential which can be used di- 
rectly to derive the shell- model effective interaction V e g. 
We have performed shell-model calculations [H C3- E3| us- 
ing this method as well as the traditional G-matrix one. 
Comparison of the results has shown that the former pro- 
vides an advantageous alternative to the latter. 

The main purpose of this paper is to give a brief sur- 
vey of the practical application of the Viow-k approach in 
realistic shell-model calculations. In section 2 we give an 
outline of the derivation of Viow-k and V g, while in sec- 
tion 3 we show a few selected results of our current study 
of nuclei around doubly magic 132 Sn, which have been ob- 
tained starting from the CD-Bonn NN potential ||. In 
section 4 we present some preliminary results of a study 
performed with different phase-shift equivalent JVJV po- 
tentials and examine their dependence on the choice of 
A. 



II. OUTLINE OF THEORETICAL 
FRAMEWORK 

As pointed out in the Introduction, we "smooth out" 
the strong repulsive core contained in the bare NN po- 
tential Vnn by constructing a low-momentum potential 
Viow-k- This is achieved by integrating out the high- 
momentum modes of Vnn down to a cutoff momentum 



A. This integration is carried out with the requirement 
that the deuteron binding energy and low-energy phase 
shifts of Vnn are preserved by Viow-k- This requirementy 
may be satisfied by the following T-matrix equivalence 
approach. We start from the half-on-shell T matrix for 
Vnn 

/>oo 

T(k',k,k 2 ) = V NN (k',k) + p / q 2 dqV N N(k',q) 

Jo 

-^T(q,k,k 2 ) , (1) 

where p denotes the principal value and fc, kl , and 
q stand for the relative momenta. The effective low- 
momentum T matrix is then defined by 



Ti ow - k (p',p,p 2 ) = Vow-k(p',p)p / q 2 dqV iow - k (p',q) 

Jo 

^-^T low ^ k (q,p,p 2 ) , (2) 

where the intermediate state momentum q is integrated 
from to the momentum space cutoff A and (p',p) < A. 
The above T matrices are required to satisfy the condi- 
tion 

T(p', Pl p 2 ) = T low ^(p',p,p 2 ); (p',p)<A. (3) 

The above equations define the effective low-momentum 
interaction Vow-k, and it has been shown |lj that they 
are satisfied by the solution: 



V\ 



OW — k 



= Q-Q' / Q+Q' 



'-Q' J QJQJQ+- 

(4) 

which is the well known Kuo-Lee-Ratcliff (KLR) folded- 
diagram expansion 0, Q, originally designed for con- 
structing shell-model effective interactions. In the above 
equation Q is an irreducible vertex function whose inter- 
mediate states are all beyond A and Q' is obtained by 
removing from Q its terms first order in the interaction 
Vnn ■ In addition to the preservation of the half-on-shell 
T matrix, which implies preservation of the phase shifts, 
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this Mow-k preserves the deuteron binding energy, since 
eigenvalues are preserved by the KLR effective interac- 
tion. For any value of A, the low-momentum effective 
interaction of equation (4) can be calculated very accu- 
rately using iteration methods. Our calculation of Vi ow -k 
is performed by employing the iteration method proposed 
in which is based on the Lee-Suzuki similarity trans- 
formation |9j. 

The Mow-k given by the T-matrix equivalence ap- 
proach mentioned above is not Hermitian. Therefore, 
an additional transformation is needed to make it Her- 
mitian. To this end, we resort to the Hermitization pro- 
cedure suggested in , which makes use of the Cholesky 
decomposition of symmetric positive definite matrices. In 
a very recent paper it has been shown that the Viow-k 
obtained through this method belongs to a family of Her- 
mitian low-momentum NN interactions which preserve 
the deuteron binding energy and are phase shift equiv- 
alent, all reproducing empirical phase shifts up to the 
cutoff momentum A. From a comparison of shell-model 
matrix elements calculated with Hermitian Vi ow -k's given 
by various methods it has turned out that they are 
approximately equivalent to each other. 

Having outlined the derivation of the low-momentum 
interaction Mow-k from Vat at, we turn now to shell- model 
calculations in which we use Vi ow -k as the input interac- 
tion. We derive the model space effective interaction Kff 
by employing a folded-diagram method, which has been 
successfully applied to many nuclei [TT| | using G-matrix 
interactions. Since Vj ow -k is already a smooth potential, 
it is no longer necessary to calculate the G matrix. We 
therefore perform shell-model calculations follo wing the 
same procedure as described, for instance, in [l2l Il3| ex- 
cept that the G matrix used there is replaced by Mow-k- 
More precisely, we first calculate the Q-box including di- 
agrams up to second order in Vi ow _k- The shell-model 
effective interaction is then obtained by summing up the 
Q-box folded diagram series using the Lee-Suzuki itera- 
tion method Q- 



III. RESULTS AND COMPARISON WITH 
EXPERIMENT 

Shell-model studies of neutron-rich nuclei in the close 
vicinity to doubly magic 132 Sn offer the opportunity to 
test the matrix elements of the effective interaction well 
away from the valley of stability. We have recently stud- 
ied |14j several of these nuclei using as input interaction 
the Viow-k derived from the CD-Bonn NN potential. We 
report here some selected results of our current study of 
nuclei with valence neutrons in the 82-126 shell, focusing 
attention on 134 Sn. This is the heaviest Sn isotope for 
which some experimental information on excited states 
is presently available. 

In our calculation we let the two valence neutrons oc- 
cupy the six single-particle levels 0h 9 / 2 , I/7/2, 1/5/21 
2p 3 / 2 , 2p 1 / 2 , and 0? 13 / 2 of the 82-126 shell. As regards 



the single-neutron energies, they have been taken from 
the experimental spectrum of 133 Sn, except that relative 
to the «i3/2 level which has not been observed. The lat- 
ter has been taken from Ref. ^jjj], where all the adopted 
values are also reported. In this connection, it should be 
mentioned that the calculations performed in this pre- 
vious work differ from the present ones mainly because 
we used there the G-matrix formalism. The results pre- 
sented in this section have all been obtained with a Vi ow -k 
derived from the CD-Bonn NN potential and confined 
within a cutoff momentum A = 2.1 fra . This value has 
been chosen according to the criterion discussed in . 
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FIG. 1: Experimental and calculated spectrum of Sn. 



From figure 1 we see that the calculated spectrum of 
13 Sn is in g ood agreement with the existing experimen- 



tal data |16l.ll7j. In figure 2 we compare with experiment 
our values of the excitation energy of the 2 + state for the 
three isotopes 128 Sn, 130 Sn and 134 Sn. The experimental 
value for 132 Sn is also reported to show the trend across 
the N = 82 magic number. The neutron single-hole en- 
ergies used in the calculations for 128 Sn and 130 Sn are the 
same as those adopted in our recent study [18j of 129 In 
and 129 Sb. 

Very recently 0, the B(E2: 0+ -> 2f ) values in 132 Sn 
and 13 Sn have been measured using Coulomb excita- 
tion of neutron-rich radioactive ion beams. We have 
calculated this B(E2) for 134 Sn with an effective neu- 
tron charge of 0.70e, according to our early study [l5| . 
The experimental and theoretical values are shown in fig- 
ure 3, where we also compare with experiment | 20j our 
B{E2;0+ -> 2f) values for 128 Sn and 130 Sn calculated 
with an effective neutron-hole charge of 0.78e |l8|. 

From figures 2 and 3 we see that our calculations repro- 
duce remarkably well the trend of both the first excited 
state and the B{E2; + — > 2+) across the N = 82 magic 
number. 
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FIG. 2: Energy of the first 2 + excited state in Sn isotopes 
around 132 Sn. The calculated results are represented by cir- 
cles and the experimental data by triangles. 
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FIG. 3: Experimental and calculated values of B(E2; + — > 
2+) for Sn isotopes around 132 Sn. 



IV. PHASE-SHIFT EQUIVALENT 
NUCLEON-NUCLEON POTENTIALS AND 
SHELL-MODEL CALCULATIONS 

The NN potential Vj ow -k is currently attracting much 
attention. On the one hand, it is being applied in various 
contexts, for instance in the study of few-body problems 
pH 122^ . On the other hand, its properties are being 
studied in detail, in particular as regards the dependence 
on the original potential model. In fact, there is evidence 
[23L l24l l25j that Viow-k's derived from different phase- 
shift equivalent potentials, such as CD-Bonn, Nijmegen 
[26| . Argonne vi$ (AV18) [27|j. are nearly identical. 

We are currently investigating 14] the dependence of 
nuclear structure results on the NN potential used to 
derive the shell-model effective interaction through the 
Vjow-k approach. Here we present some preliminary re- 
sults obtained for the nucleus with two valence protons 
134 Te, which provides a good testing ground for this in- 
vestigation. We consider the three phase-shift equivalent 
NN potentials Nijmegenll, AV18 and CD-Bonn. As re- 
gards the cutoff A, we let it vary from 1.5 to 2.5 fm -1 . 



TABLE I: Energy levels (in MeV) of 134 Te . Predictions by 
different NN potentials are compared with experiment. 





Nijmll 


AV18 


CD-Bonn 


Expt 


E ga 


-20.793 


-20.778 


-20.778 


-20.560 


E(2+) 


1.342 


1.325 


1.315 


1.279 


E(A+) 


1.611 


1.599 


1.603 


1.576 


E(6+) 


1.741 


1.726 


1.742 


1.691 



In all cases we compare the calculated spectrum of 134 Te 
with the experimental one up to about 4.5 MeV exci- 
tation energy (this includes sixsteen excited states) and 
calculate the rms deviation eryfj. 
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FIG. 4: Behavior of the rms deviation a relative to the spec- 
trum of 134 Te as a function of A for different NN potentials. 
See text for details. 

In figure 4 we show the behavior of a as a function of A 
for the three potentials. We see that the curves relative to 
Nijmegenll (Nijmll) and AV18 practically overlap each 
other while that for the CD-Bonn potential has a rather 
different pattern. In particular, the minimum of a for 
CD-Bonn is at A ~ 1.7 fm -1 while for the other two 
potentials it lies at A ~ 2 fm -1 . The minimum value of a 
for the three potentials is however almost equal and also 
the energies of the various states are practically the same. 
By way of illustration, we report in table 1 the ground- 
state energy of 134 Te (relative to doubly magic 132 Sn) and 
the excitation energies of the first three positive-parity 
yrast states. 

From figure 4 it also appears that all three curves are 
rather flat around the minimum. As a consequence, the 
quality of agreement between theory and experiment does 
not change significantly for moderate changes in the value 
of A. More precisely, for the CD-Bonn potential a re- 
mains below 150 keV for values of A between 1.5 and 1.9 
fm -1 while for the other two potentials this occurs for A 
between 1.8 and 2.2 fm -1 . 

The main conclusion of this preliminary study is that, 
allowing for limited changes in the value of A, nuclear 
structure results obtained from Vi ow -k's extracted from 
different NN potentials are practically independent of 
the input potential. It therefore appears that the low- 
momentum interaction Vj ow -k gives an approximately 
unique representation of the NN potential. This is quite 
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in line with the conclusions of |28j . 

At this point, a comment on the results presented in 
this section is in order. As mentioned in section 2, we 
include in the Q-box all diagrams up to second order in 
Mow-k- In the calculation of these diagrams we have in- 
serted intermediate states composed of particle and hole 
states restricted to the two major shells above and be- 
low the Fermi surface. However, as a development of our 
study, we are currently investigating the effect of increas- 
ing the number of intermediate states. The first results 
indicate that the minimum value of a occurs at values of 



A which are somewhat larger than those reported here. 
Our final findings will be the subject of a forthcoming 
publication. 
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We define a = {(l/N d ) [E exp (i) - 
Nd is the number of data points. 



E calc (i)] } ' , where 



